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Abstract 

Three exactly solvable Hamiltonians of complex structure, are studied in 
the framework of a semi-classical approach. The quantized trajectories for 
intrinsic coordinates correspond to energies which may be classified in collec- 
tive bands. For two of the chosen Hamiltonians the symmetry SU(2)®SU(2) 
is the appropriate one to classify the eigenvalues in the laboratory frame. 
Connections of results presented here with the molecular spectrum and 
Moszkowski model are pointed out. The present approach suggests that 
the intrinsic states, which in standard formalisms are heading rotational 
bands, are forming themselves "rotational" bands, the rotations being per- 
formed in a fictitious boson space. 

I. INTRODUCTION 

Collective motion in nuclei has constantly been an exciting subject 
for both experimentalists and theoreticians. Some phenomenological mod- 
els have been successfully in interpreting the data for energies and transition 
probabilities for certain ranges of atomic mass. The long series of publi- 
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cations has been opened by the hquid drop model [T| which was later on 
improved by including anharmonicities, deformed equilibrium shapes or by 
extending the symmetries for the intrinsic motion ^|-^] . Some phenomeno- 
logical models interpret the properties of collective states of positive parity 
in terms of quadrupole degrees of freedom. Interacting boson model (IBA) 
points out the fact that the monopole degrees of freedom should be 
included which results in obtaining a conceptually new formalism for the 
description of nuclear dynamics. All phenomenological models define no- 
tions like rotational bands, equilibrium shapes, nuclear phases based on 
the behavior of the system in the intrinsic frame. The coherent state model 
(CSM) exploits the classical features of collective states of high angular mo- 
mentum and provides a realistic description for deformed and transitional 
region in terms of angular momentum projected states from coherent states 

i- 

Some microscopic theories were developed by paralleling the aims 
of the phenomenological models. Indeed, the nuclear properties are de- 
scribed in terms of individual degrees of freedom. In this sense to define 
the optimal collective coordinates for a many body system, was always a 
central field of activity. Microscopic formalisms use approximations like 
Bogoliubov-Hartree-Fok, Random Phase Approximation (RPA), boson ex- 



pansion technic, higher order RPA, extended shell model, etc. [10|, which 
need to be tested and define the circumstances under which they work or 
not. For such purposes very often one uses completely solvable models 
where the results of various approaches are compared with the exact re- 
sult. Most known schematic and solvable models are those of Moszkowski 
2|, Lipkin-Meshkov and One Level Pairing Hamiltonian |12,13|. Also 



phenomenological models proposed by Wilets and Jean [^, Davydov and 
Filippov are even nowadays used to get a reference framework which 
allows for an interpretation of the results obtained with more sophisticated 



models. 

In this paper we address the question whether the motion of intrinsic 
collective coordinates can be described by the irreducible representation of a 
symmetry group. The chosen model Hamiltonians have a complex structure 
being of fourth order in quadrupole bosons. They have however a common 
feature that all of them are, in a semi-classical picture, completely solvable 
and therefore analytical expressions for energies are obtained. 

The model Hamiltonians are explicitly treated according to the fol- 
lowing plan: In Section 2 a fourth order Hamiltonoian consisting in a 
quadratic quadrupole boson term plus a fourth order term which is the 
square of the second order boson invariant. Two sets of parameters are con- 
sidered, one yielding, after dequantization, an effective potential with one 
minimum and an other one which produces a pocket potential. Solutions 
are presented for both cases. Energies were organized in bands like standard 
bands are in the laboratory frame. The difference is that here the angular 
momentum is defined by an TZ^ group acting in a fictitious space. In Section 
3 the two body boson interaction is of a multipole-multipole type with the 
multipolarity 2^,k = 0,2,4. In the intrinsic frame the classical Hamilto- 
nian has energies which after quantization might be classified by L,M and 
therefore the corresponding states are SU(2) states. This pictures implies 
that in the laboratory frame the symmetry SU(2)(8)SU(2) might be used. 
In Section 4, the two body boson interaction is of a quadrupole-quadrupole 
type. Neglecting the correlations with the excited collective bands the clas- 
sical Hamiltonian is identical with those proposed by Moszkowski |22] and 
therefore a simple expression for the corresponding spectrum is possible. 
Final conclusions are drawn in Section 5. 



3 



II. THE FIRST EXAMPLE OF AN EXACTLY SOLVABLE HAMILTONIAN 



In this Section we analyze the classical behavior of a fourth order 
boson Hamiltonian with a particular structure. The third order term is 
missing and the fourth order one depends exclusively on quadrupole coor- 
dinates but on their conjugate momenta: 

Hi = ^1(6^62)0 + A2 [(6^6^)0 + (6262)0] + ^4P'. (2.1) 

Here h^^ (^2/*) denotes the creation (annihilation) quadrupole boson opera- 
tor while P is the second order boson invariant: 



We consider the variational equation 

6 j\ij\Hi-ih-^\i;)dt' = (2.3) 

where | denotes the following coherent state 

I ^) = exp[zobto - z^b2o + Z2{bt2 + bt_2) - 4{b22 + 62-2)] I 0) (2.4) 

The vacuum state for the quadrupole boson operators 62,^ is denoted by 
I 0); zo and Z2 are complex functions of time and play the role of classical 
coordinates. To simplify the notations, hereafter the quadrupole boson 
operators will be denoted by 6^ omitting the index specifying the angular 
momentum carried by the boson operators. 
The coordinate transformation: 

fc+2 

qi = 2 i Re{Zk), 

Pi = 2^Im{zk), k = 0,2, ^ = ^, (2.5) 

brings the classical equation of motion to a canonical form: 

d Hi 



dUi 



dpk 
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-Pk, 

Qk- (2.6) 



Here TCi stands for the average value of the chosen Hamiltonian on the 
coherent state | while "dot" denotes the time derivative operation. 

Hi ^ (V I ^1 I ^) = y (p? + pI) + ^iql + ql) + D{ql + qlf. (2.7) 

The coefficients involved in the expression of 7ii are related to those defining 
the model Hamiltonian by the following equations: 

A=^{Ai+2A2) + UAi, 
^' = ^(^1-2^2), 

D = 4.A^. (2.8) 

For what follow, it is useful to introduce the polar coordinates: 

(71 = r cos 6* , q2 = rs\n9. (2-9) 

In terms of the new coordinates the classical energy function has the ex- 
pression: 

m = ^(r^ + rH^) + ^r^ + ^r\ (2.10) 

From the above equations we see that the classical system associated to the 
boson Hamiltonian is a particle which moves in a plane due to the force 
determined by a potential invariant with respect to rotations around an 
axes perpendicular to the motion plane: 

V{r) = ^r' + ^r\ (2.11) 

The potential energy is plotted in Fig.l for three sets of parameters (A,D). 
The parameters were chosen so that three distinct situations are pointed 
out, namely when the equilibrium shape is spherical, deformed and meta- 
stable. Using the equations of motion ( p.6| ) one can prove that 

£3 = , Wi = 0. (2.12) 

where £3 is defined by the following expression: 



-Cs = ■7^{qiP2 - q2Pi) = 



Let us define another two classical functions on the phase space, spanned 
by the coordinates {qi,Pi,q2,P2) 

Ci = \{{ql + pi-ql-pl) 

^2 = 2 +P1P2)- 

Given two function /i and /2, defined on the phase space, their Poisson 
bracket is defined as: 



{/l,/2}= E 



fc=l,2 



dhdf2 dhd f2 



d Qkd pk d pkd qk 



The classical functions Ck obey the following equations: 

{£i,£2} = fiCs, 
{£2, A} = ^^1, 
{£3,£i} = hC2. 

In virtue of Eq.( |2.16D the set of functions Ck with the Poisson brackets as 
multiplication operation, form a classical SUc{2) algebra. Moreover they 
could be obtained by averaging on | the generators of a boson SUb{2) 
algebra defined with the boson operators ftg", b^2- 

Ck = {ip \ Lk \ ip)]k = 1,2,3, 



Li = - 



26l,6o-(4 + ^-2)(^2 + fe- 



L3 



h 



2V2 
h 

2y/2i 



bUb2 + b-2) + {bl + bl2)bo 



bl{b2 + b^2) - {bl + bU)bo 



The first equation ( 2.17] ) and the correspondence between commutators and 
Poisson brackets 
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define a homeomorphism of tlie boson and classical algebras generated by 
{Lk}k=i^2,3 and {/^fc}fc=i,2,3 respectively. Note that the boson SUh{2) alge- 
bra does not describe the rotations in the real configuration space but in a 
fictitious space. The conservation law expressed by ( p. 12 ) is determined by 



the invariance against rotation around the 3-rd axis in the fictitious space 
mentioned above: 

[H,,Ls]=0. (2.19) 

Since the classical system is characterized by two degrees of freedom and, 
on the other hand, there are two constants of motion 

ni = E , C3 = ^, (2.20) 

the equations of motion are exactly solvable. Indeed, inserting the constants 
of motion in Eq. ( 2.10| ) the resulting differential equation 



^ + yeff{L;r) = E, (2.21) 

with 

^e//(L;r) = ^ + f^ + fr^ (2.22) 
provides the time variable as function of x = by 

^ 1'^ (2.23) 



D y ~^^^y ^ D 

Note that L has the meaning of the third component of the angular mo- 
mentum in the space spanned by the coordinates qi,q2,q3 with a third 
coordinate which might be associated in a more complex situation to an 
additional degree of freedom. 

In our numerical application we considered two sets of parameters: 

I) A' = OMMeV, A = 3MeV, D = OAMeV, 

II) A' = O.OlMeV, A = 3MeV, D = -OMMeV. (2.24) 



The two situations are represented pictorially in Figs. 2,3 for several val- 
ues of L. One notes that while the first effective potential has only one 
extreme value, a minimum, which depends on angular momentum, in the 
case labeled with II, there are two extremes, one minimum and one maxi- 
mum. While the minima are depending on angular momenta, maxima are 
almost independent. Moreover there is a critical value for angular momen- 
tum where the two extremes get unified into an inflexion point. Below the 
critical angular momentum the effective potential has a pocket shape and is 
similar to that obtained in the study of heavy ion collision with two centers 
harmonic potential. Equation ( 2.21| ) suggests that the motion is allowed 



only for energies obeying the restriction: 

E > Vjff (2.25) 
For a given pair of (E,L) the intervals of r were the motion takes place are: 
r-i < r < ra, and r > rg for E < V^'''', 

r>ri for E > V^^^, (2.26) 

were denotes the values of r were E = V^jf. The minimum and maximum 
values of effective potential are denoted by VJpJ" and ^Jff^j respectively. 
Classical trajectories may evolve, depending on the initial conditions, either 
on a finite trajectory in the interval ri < r < r2 or on an unbound one with 
r > r-^. As we shall see, in the first interval the classical motion could be 
quantized while in the second one the system undergoes a fission process. 
In our semi-classical quantization procedure the system cannot go, through 
a tunneling effect, to the unstable region. However, in a full quantum 
mechanical description the wave function describing the system inside the 
pocket region, is spread also to the region r > r^. 

Since t is a real quantity the integration interval in (2.23) is chosen 



so that inside it the polynomial is positive. As shown in Figs. 4 and 5 these 
intervals are: 



e-2 < X < e-i for case I, ei < x < 62 for case II. 

This suggests that a possible initial condition is: 
' 62, case I 

xq = x(0) = < ei. case II, bound trajectory 

. 63, case II, unbound trajectory . 
The integral involved in Eq. ( 2.23| ) can be analytically performed [16| 



and the final result for time is 



^\A'D\{e^-ei\ 



sn ^{u,k), 



(63 - ei){x - 62) 63-62 

-7 r , k = \ , tor the case i 

(63 - e2)(x - ei) V63-61 



X — 


61 


62 - 


61 


63 - 


61 


X — 


61 



k = W , for the case II, 61 < 2; < 62 



62 


- 61 


63 


- 61 


62 


- 61 


63 


- 61 



k = J , for the case II, x > 63. 

V 63 - ei 

Here sn{u,k) denotes the elliptic sinus function. Its argument "u" depends 
on 'x', the upper limit of the integral defining the time as well as on the 
coordinates 61, 62, 63 where the polynomial lying under the square root sym- 
bol vanishes. Actually these are obtained by squaring ri, r2, r3 respectively, 
the turning points of the effective potential mentioned above. The roots 6^ 
are ordered as follows: 

ei < 62 < 63. 
This function can be inverted 

and finally the coordinate r is expressed as a function of time. The explicit 
expression for the coordinate x as a function of time is given in Appendix A. 
Using this result in connection with the second equation (2.13), the equation 
for the coordinate is solved 

A'L 



n^git') 
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-dt' = h{t) 



and thereby the classical trajectory is fully determined. 

An alternative way to derive the implicit equation of the classical 
trajectory is described in Appendix A. 

For the variable x, the motion is periodical with the period: 

where 2F1 denotes the confluent hyper-geometric function with the argu- 
ment /c^, where k is defined by Eq.( |2.29| ). 

Note that the period depends on the energy E, by means of the 
roots ei , 62 , 63 . Energy is discretitized through the quantization equation 
of Bohr-Sommerfeld type: 



The function appearing in the left hand side of the above equations is plotted 
in Figs. 6, 7 as function of E for the two particular potentials considered 
here. One remarks the fact that this integral depends almost linearly on 
energy, irrespective of L. However the slope of these straight lines depends 
on the potential characterizing the motion of the system. 

Quantized energies for n=0,2 and L=even and those with n=l and 
L > 2 obtained for the pocket like potential are plotted in Fig. 8. One 
should note that although these states describe the intrinsic degrees of free- 
dom, they might be organized in bands as it happens with the states char- 
acterizing the whole set of degrees of freedom. In the case II the bands 
are finite since trajectories with energy larger than VJ^j^ cannot be quan- 
tized, describing unbound systems. Obviously the number of states which 
could appear inside the pocket depends on angular momentum. The larger 
angular momentum, the smaller the number of the bound states. We men- 
tion again that in a pure quantum mechanical treatment these states are 
quasi-bound. In contrast with this situation, for the case I, the bands are 
all infinite. 




(2.34) 
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Prom Figs. 6-8, it results that for a given value of angular momentum 
the energy spacings of states with different values for the quantum number 
n is almost constant. This feature is determined by the relative values of the 
coefficients defining the model Hamiltonian. For this particular situation, 
the harmonic approximation seems to be quite reasonable. Indeed, expand- 
ing the effective potential around its minimum and ignoring the terms of 
order higher than two, one finds for the harmonic frequency the value: 



(2.35) 



Thus, the lowest three bands can be defined in an approximative way by: 

4'^=K7f(^,^o),i> = 0,2,4,6,... 
4''^ = ^eF(^'^o)+u;L,i: = 2,3,4,... 

= V,JfiL,ro) + 2uL,L = 0,2,4,6,... (2.36) 

Exact energies of the levels in the three bands are given in Fig. 8 as 
function of angular momentum L for the case II. From Fig. 7 one may see 
that the number of the bound states in the effective potential 

characterizing the case II depends on the value of the angular momen- 
tum. This dependence is seen more clearly in Fig. 9. From there one sees 
that the number of bound states varies from 18, for L = 0, to 0, for L > 46. 
Concluding this section, although the starting Hamiltonian is anharmonic 
and has a complex structure we derived analytical solutions for the classi- 
cal trajectories. The corresponding energies where quantized using a Bohr 
Sommerfeld type quantization condition. Energy levels are labeled by two 
indices, n,L, one provided by the quantization procedure and the other one 
by a conservation law. We fixed, conventionally, the constant value of £3 
to be equal to L/2. However the classical trajectory is not invariant to 
the fictitious group SUc{2). Therefore it is deformed being a mixture of 
components of different L. This can be easily seen if we evaluate JC^ for a 
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trajectory of energy E and projection of angular momentum on z-axis equal 
to M: 



+ 



4:A 



/2 



E + {A'- Ay - jr^ 
E + {A - A)r^ - 



n 2 



2li 



2 • 



From this equation it is manifest that is not a constant of motion since 
this quantity depends on r and therefore changes its value when the system 
moves from one point to another. 

The problem treated in this Section is a real challenge for looking for 
a quadrupole boson Hamiltonian which might yield in the intrinsic frame a 
motion which admits both C^. and as constants of motion. This would 
allow us to classify the intrinsic motion with a SU{2) symmetry. 

III. ANOTHER SOLVABLE BOSON HAMILTONIAN 



Here we study the fourth order boson Hamiltonian: 



M J=0,2,4 



Along the time , this Hamiltonian has been used by several authors to de- 
scribe the rotational ground band. First authors were Das, Dreizler and A. 



Klein |17| who treated this Hamiltonian in a particular basis and obtained 
an analytical formula for the angular momentum dependence of the yrast 
states energies. Later on this boson number conserving Hamiltonian was 



used by lachello Ref. |18,19| to formulate the first version of the interacting 
boson approximation. In |2^, this Hamiltonian was averaged on a angu- 
lar momentum projected state obtained from an axially deformed coherent 
boson state to approximate the energies from a rotational ground band . 
In the vibrational limit the empirical formula of Ejiri |25] was rigorously 
derived. 



12 



As shown in Appendix B, the Hamiltonian H2 can be written in an 
equivalent otherwise more convenient form: 

H2 = {l + e)N2 + PNI + 5JI + aib^b^)oibb)o 

with the coefficients defined there. The notation N2 is used for the 
quadrupole boson number operator while Jf stands for the total angular 
momentum squared carried by the quadrupole bosons. The classical energy 
function associated to H2 has the expression: 

n2 = A{ul + V^ + 2u\ + vl)) + B{ul + + 2u\ + vl)f + C{uqV2 - U2VQf , 

with the coefficients A, B, C given in Appendix B. 

The classical equations of motion are obtained from the time de- 
pendent variational equation ( |2.3D by replacing Hi with H2. They have 
a canonical form with respect to the coordinates {qkiPk}k=i,2 defined by 
Eq. ( |2.5| ). In terms of canonical coordinates the energy function has the 
expression: 

^2 = ^{ql + ql+pl+pl) + j{ql + ql + pI + P^f + j{qiP2 - q2Pif ■ 

Since the boson Hamiltonian commutes with L3 and there are two in- 
dependent constants of motion: 

AA2 = (^liVal^), £3 = (^|4|*) 
This results from the equations: 

+ = 0, 

d 

— (giP2 -P192) = 0, 

implied by the equations of motion. Using these results one finds out that 
7i2 is a constant of motion. On the other hand, considering the expressions 
of the SUc{2) generators one finds that: 
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\{<ll + <i+pl+pl^ . (3.7) 

and according to the previous equation this is a constant of motion. The 
classical energy can be quantized in two equivalent ways, namely either one 
quantizes the angular momentum by fixing the constants of motion such 
that: 

= h^L{L + 1), £3 = hM, with -L<M <L, L,M integers. (3.8) 

or by quantizing the classical action 

J (QiPi + Q2P2)dQdp = 2Trfin, n — positive integer (3-9) 

and the third component of angular momentum. For the first option the 
result for the quantized energy is: 

Elm = 2A^L{L + 1) + 4BL{L + 1) + ^M^, (3.10) 

while in the second case the result is: 

E„,M = A{n + l)+B{n + lf + ^M^ (3.11) 

In this equation we introduced the zero point motion for the plane oscillator 
quanta although the semi-classical quantization in not able to account for 
it. To compare the two alternative expressions for the quantized energy it 
is convenient to make the following approximations, which work quite well 
for large quantum numbers: 

(n+l)2«n(n + 2), ^L(L+1)«L + ^ (3.12) 

In this way the two expressions for energy are identical provided ^ = f ■ It is 
worth mentioning that the Hamiltonian considered in this section was used 
in Ref. 17 in the boson basis {\NvaJM)}, where the quantum numbers 
involved are the number of bosons, seniority, missing quantum number , 
angular momentum and its projection on z-axis, with the restriction N = 
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V = ^ for the ground band states. By contrast, the connection between 
the new angular momentum quantum number L and the number of quanta 
in plane is L = ^. This symmetry obeyed by the plane oscillator was 



exploited by Moszkowski in Ref. |22| where a schematic solvable many body 
Hamiltonian was found which describes fairly well the main features of the 
competition between individual and collective degrees of freedom. 



The equation ( ^.10 ) shows that for the model Hamiltonian considered 



in this section the energy of the intrinsic degrees of freedom can be classified 
by the quantum numbers (L,M) and correspondingly the quantized states 
by the irreducible representations of the fictitious SU},{2) group. Since 
energies depend on M, the states are deformed despite the fact L is a good 
quantum number. From the classical energy function it is clear that one 
deals with a plane oscillator. For even values of "n" one may define a one 
dimensional oscillator whose number of quanta is determined as 2u = n. 
Also we fix the projection of angular momentum on z-axis by the condition 



£3 = y/L{L + 1). Under these circumstances the quantized energy gets the 
expression: 

E,,L = 2A{u + i) + iB{i. + i)2 + ^L(L + 1). (3.13) 

Apart from an additive constant term this equation coincides with those 
used by Erb and Bromley to fit the spectrum of ^^C C system [^,21]. 



Using a more complex structure for the starting model Hamiltonian, 
a coupling of rotational and vibrational degrees of freedom is possible. Such 
a coupling has been described by lachello p7|-p9t| in the framework of an 
algebraic model. 

Equations derived in this Section for energy refer to the intrinsic de- 
grees of freedom. Passing to the laboratory frame, the total energy accounts 
also for the rotational degrees of freedom. Assuming that the intrinsic de- 
grees of freedom and Eulerian angles, defining the position of the intrinsic 
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frame, are only weakly coupled, the total energy can be approximated by: 

C 



EjLM = 2A^L{L + 1) + ABL{L + 1) + ^M^ + 5J{J + 1). (3.14) 

Concluding this Section one may say that the energies of the exactly 
solvable Hamiltonian can be classified by a SU(2)(g)SU(2) symmetry, one 
factor describing the motion of intrinsic degrees of freedom while the other 
one taking care for the motion of Eulerian angles. 

In the next section we consider a boson Hamiltonian to which it 
corresponds a classical energy function which is similar to the mean field 
underlying the model proposed by Moszkowski [p^] . 

IV. A QUADRUPOLE-QUADRUPOLE BOSON HAMILTONIAN. 

In this section we study a quadrupole boson Hamiltonian which, in the 
intrinsic frame is closely related to the schematic Hamiltonian introduced 
by Moszkowski p^, long time ago. This Hamiltonian has been used by 



several authors to test various many-body approaches |23,24|. For the sake 
of completeness we present first the main ideas underlying the Moszkowski 
model (MM). 



A. Brief review of Moszkowski model 

The MM model considers a system of nucleons moving in a mean 
field, consisting in a two dimensional oscillator potential and a spin-orbit 
term, and interacting among themselves through a quadrupole-quadrupole 
interaction. Let us consider first the one body Hamiltonian: 

^ Hho + Hso = ^{pI + pD + ^{x^ + y^) - Cls. (4.1) 
It can be checked that Hy^o commutes with the quasi-spin operators: 
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1 

4 

_ 1 

ty-2 

tz = ^Iz = ^{xPy-yPx)- 



Muoxy + 
1 



Mloq 



PxPy 



and therefore its eigenstates can be classified by the irreducible representa- 
tion of the SU{2) group generated by {tA;}A;=a;,j/,2- 

P\Nma) = —{— + l)\Nma) = t{t + l)\Nma), 
tz\Nma) = —m\Nma). 

Here N denotes the total number of quanta and m the azimuthal quantum 
number. The spin is perpendicular on the (x,y) plane, the two possible 
orientations being specified by a{= it). It is worthwhile to notice that in 
terms of stretched coordinates: 

y 

\zl 



y' 

\z'J 



iMujQ 

n 



the components of quasi-spin operator are formally identical with those of 
angular momentum Ck, defined by Eqs. (2.13), (2.14). However the two 
sets of operators act on different spaces. 

Let us consider a many body system moving in the mean field de- 
scribed by Hgp and interacting by the QQ force, which in plane acquires a 
very simple form: 

1 , . X 



QQ 



where T denotes the operator acting on the many body states: 



{Nma\t^\Nma)cl 



m,m' 



One distinguishes two limiting cases. When C=0 one obtains the two di- 
mensional version of the Elliott model [EH] which is suitable for describing 
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the collective rotations and quadrupole vibrations of a many body system. 
The other regime, when the long-range interaction is missing (X=0), sim- 
ulates the shell model description in realistic situations. The intermediate 
situations can be covered by a smooth variation of the two strength param- 
eters C and X. Renormalizing the mean field due to the two body interac- 
tion one obtains a single particle Hamiltonian which in terms of stretched 
coordinates looks as: 

H'sv = |(^" + y" +p.''+Py'') - + y" +p.'' +Py''f + §{^'Py' " y'p.'f- (4.7) 

As we shall see, this expression for the renormalized mean field is very 
useful. 



B. Two body interaction for quadrupole bosons. 

In this subsection we shall treat semi-classically the following 
quadrupole boson Hamiltonian: 

H3 = nN2-^Yl Q2^Q2-M(-)^ (4-8) 
where denotes the quadrupole operators defined as: 

Q2^ = VQ{b^b)^^. (4.9) 

Averaging this boson Hamiltonian on the coherent state ( [2.4] ) and approx- 
imating the average for the two body term by J2^{Q2fi){Q2-fi{—)^), one 
obtains the following expression for the classical energy function: 

Comparing the energy function with the mean field corresponding to 
Moszkowski model we realize that they are identical if we accept the follow- 
ing relationships between the strengths involved in the two Hamiltonians: 

n = uj, -F = —X. (4.11) 
' 7 32 ^ ' 
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Therefore the spectrum of the quantized intrinsic Hamiltonian is: 

Elm = 2n^L{L + 1) - yF[L(L + 1) - M^]. (4.12) 

The competition between the individual and collective feature in 
Moszkowski model corresponds to the interplay of harmonic and anhar- 
monic terms in the boson interacting model. 

Again we pointed out the possibility to classify the states describing 
the intrinsic degrees of freedom by the irreducible representation of an SU(2) 
group. 

V. CONCLUSIONS. 

Three boson Hamiltonians of complex structure otherwise exactly 
solvable, were semi-classically treated. The first is a boson number non- 
conserving Hamiltonian. Its classical trajectories are analytically expressed 
in terms of elliptic functions. Closed orbits are quantized and explicit ex- 
pressions for discrete energies are obtained. They are organized in rota- 
tional bands following the traditional scheme of the liquid drop model. The 
difference is that here the states describe the motion of intrinsic degrees 
of freedom and moreover the angular momentum components generate a 
rotation group acting in a fictitious space. 

The second Hamiltonian commutes with the quadrupole boson num- 
ber operator and yields in a classical framework a spectrum which is classi- 
fied by angular momentum and its projection on z-axes. In the laboratory 
frame the spectrum is classified by the symmetry SU(2)(8)SU(2). The simi- 
larity with the molecular spectrum used for studying the quasi-bound states 
in ^^C +^'^ C system was pointed out. 

The third Hamiltonian includes, as a two-body boson interaction, 
the separable quadrupole-quadrupole interaction. The classical image of 
this Hamiltonian is similar to the schematic Hamiltonian introduced by 
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Moszkowski to study the interplay of individual and collective degrees 
of freedom. Its spectrum is also classified according to the symmetry 
SU(2)®SU(2). 

Concluding, the present paper points out a certain SU(2) symmetry 
which allows for a complete description of the intrinsic dynamic variables, 
13 and 7. In a subsequent publication we shall investigate the question 
whether this symmetry is observed experimentally or not. Our optimism is 
supported by the fact that the considered Hamiltonians were already used 
in connection with different approaches to describe some of available data. 



VI. APPENDIX A 



From Eq ( 2.291 ) one can obtain analytical expressions for the coordi- 
nate X as function of time. 

When the effective potential has only one extremum, the case I, the 
result is: 



x{t) = ei + 



(ea - ei)(e2 - ei) 



(es - ei) - (es - 62)571,2 



\A'D\ 



As for the potential exhibiting a pocket structure, one obtains: 
x{t) = ei + (62 - 6i)sn^ 



J^(6,-60f-/^^ 



2h^ 



63 - 61 



, ei < x < 62, 



sn 



x{t) = 63 + (63 - 62)- 



^(e3-ei)t;yiS. 

Now we describe briefly how to obtain the implicit equation for the 
classical trajectory. Inserting the expression of the time derivatives given 
by Eqs. and (|2]2^) in the identity: 



, X < 63. 



M _ dB_dr_ 
dt dr dt 



one obtains: 



(A.l) 



(A.2) 



(A.3) 
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-2Dh^ 



dy 



The integral involved in the above equation can be analytically performed 

and the final result is: 

■^n (^arcsin{u), k"^^, kj + -^^^sn~^[u, k), case I, 62 < x < 63 



X = r . 



(A.4) 



L 



2\A' 



hV {e3-ei)\D\ 



n larcsin{u), kj , case II, ei < x < 62 



Ularcsin{u),-^^,k] + 



e3V^^=iW£^ g„-i( fe) case ll,x>e3) 



(A.5) 



The arguments u and k depend on the interval and are those defined by 
Eq. (2.29). We denoted by 11 the elliptic function of third rank defined as: 

^{(p, n, k) 



(1 + nsi'n?a)y/l — k'^sin^a 



(A.6) 



VII. APPENDIX B 



The anharmonic terms involved in Eq(^]^) defining the Hamiltonian 
H2 can be easily expressed in terms of boson number and total angular 
momentum operators: 



^ = -^\^6^6t)o(66)o + ^VbN2{N2 - 1) - ^^('^2 - QN2) 

^ = ^(6^6^)0(66)0 + \n2{N2 - 1) + ^(J2 - 6iV2). 



(6t6t)2(66)2 
(6t6t)4(66)4 

Using these expressions in connection with Eq. (|3.l| ) one obtains H2 in the 
form given by (|3.2| ) with: 



a 



Co - '^V5C2 + ^C4, 



/3 = -iV5C2 + ^C4, 
35 7 

^ = ^V5C2 - |C4, 

35 7 





Averaging H2 on the coherent state (2^) one obtains the classical energy 
from Eq. ( |3.4[) where the following notations have been used: 



(B.l) 



(B.2) 



21 



A = j + e + P + 6S, 
B = + 

C=-f. (B.3) 



22 



REFERENCES 

[1] A. Bohr, Mat. Fys. Medd. Dan. Vid. Selsk. 26 (14) (1952); A. Bohr and B. Mottelson, Mat. 
Fys. Medd. Dan. Vid. Selsk. 27 (16) (1953). 

[2] A. Faessler and W. Greiner, Z. Phys. 168 (1962) 425; 170 (1962) 105; 177 (1964) 190; 
Faessler, Nucl. Phys. 59 (1964) 177; A. Faessler, W. Greiner and R. K. Sheline, Nucl. Phys. 
62 (1965) 241; 70 (1965) 33. 

[3] G. Gneuss, U. Mosel and W. Greiner, Phys. Lett. B 30 (1969) 397. 

[4] M. Seiwert, P. Hess, J. Maruhn and W. Greiner, Phys. Rev. C23 (1981)2335. 

[5] W. Greiner and J. M. Eisenberg, Nuclear Theory, vol. I, Nuclear Models. 

[6] S. G. Lie and G. Holzwarth, Phys. Rev. C 12 (1975) 1035. 

[7] P. O. Lipas, T. Honkaranta, Phys. Scripta 13 (1976) 339. 

[8] A. A. Raduta, V. Ceausescu, A. Gheorghe and R. M. Dreizler, Phys. Lett. B 99 (1981) 
444; Nucl. Phys. A 381 (1982) 253. 

[9] A. Arima and F. lachello, Ann. Phys. 99 (1976) 253; 111 (1978)201; 115 (1978) 145. 

[10] See for example P. Ring and P. Shuck, The Nuclear Many-Body Problem, Springer, 1980. 

[11] H. J. Lipkin, N. Meshkov and A. J. Glick, Nucl. Phys. 62 (1965), 188, 199, 211 

[12] A. K. Kerman, R. D. Lawson and M. H. Macfarlane, Phys. Rev. 124 (1961) 162. 

[13] R. K. Richardson, Phys. Rev. 141 (1966) 949. 

[14] L. Wilets and M. Jean, Phys. Rev. 102 (1956) 788. 

[15] A. S. Davydov and G. F. Filippov, Nucl. Phys. 8 (1958) 788. 

[16] L S. Gradshtein and L M. Ryzhik, Table of Integrals, Series and Products, Academic Press, 
1965, New York, pages 904,917, 919. 

[17] T. K. Das, R. M. Dreizler and A. Klein, Phys. Rev. C2 (1970) 632. 

23 



[18] F. lachello in Int. Conf. on Nuclear Structure and Spectroscopy, Amsterdam, 1974. 

[19] A. Arima and F. lachello, Phys. Rev. Lett. 35 (1975) 1069, Ann. Phys. (N.Y.)99 (1976) 
253. 

[20] A. A. Raduta and R. M. Dreizler, Nucl. Phys. A 258 (1976) 109. 

[21] See for example Walter Greiner, Jae Y Park and Werner Scheid, Nuclear molecules. World 
Scientific, 1995, page 104. 

[22] S. A. Moszkowski, Phys. Rev. 110 (1958) 403. 

[23] A. A. Raduta, Phys. Rev. C 51 (1995) 2973. 

[24] E. R. Marshalek, Ann. Phys. (N.Y.)143(1982) 191. 

[25] H. Ejiri, M. Ishihara, M. Sakai, K. Katori and T. Imamura, J. Phys. Soc. Jap. 24(1968) 
1189. 

[26] K. A. Erb and D. A. Bromley, Phys. Rev. C23 (1981) 2781. 
[27] F. lachello. Chemical Physics Letters,78 (1981) 581. 

[28] O. S. Van Roosmalen, A. E. I. Dieperink and F. lachello. Chemical Physics Letters 85 
(1982) 32. 

[29] F. lachello and R. D. Levine, J. Chem. Phys. 77 (1982) 3046. 
[30] J. P. Elliott, Proc. Royal Soc. London Ser. A 245 (1958) 128. 



24 



FIGURES 




FIG. 1. The potential energy involved in Eq. (2.11) is plotted as function of r for three sets 
of parameters (A,D). 
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FIG. 3. The same as in Fig. 2, but for the case II 
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FIG. 4. The third order polynomial standing under the square root symbol in Eq. (2.23), 
multiplied with sign{—A'D) is plotted as function of x = r^, in the interval x > 62 where 
ei < 62 < 63 denotes its roots. The parameters correspond to the case I, listed in Eq. (2.24). 
The plotted function is positive in the interval 62 < x < 63. 
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FIG. 5. The same as in Fig. 4 but for the case II. The plotted function is positive in the 
interval ei < x < ei 
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FIG. 6. The function 1{E), defined by Eq. (2.34) is plotted as function of E for several 
angular momenta. The parameters involved are those specified in Eq. (2.24) by the label I. 
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E[MeV] 



FIG. 7. The same as in Fig. 6, but for the case II 
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FIG. 8. The energies associated to the motion of intrinsic degrees of freedom are classified 
in rotational bands in the following way. The lowest band is characterized by n=0 and even 
values for L. The second band is similar to the gamma band and corresponds to n=l and L > 2. 
The third band is similar to the band /3 and has n=2 and even angular momenta. 
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FIG. 9. The number of bound states in the pocket like potential is plotted as function of 
angular momentum. 
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